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ABSTRACT 

We present a high-resolution simulation of globular cluster formation in a galaxy merger. For the 
first time in such a simulation, indiyidual star clusters are directly identified and followed on their 
orbits. We quantitatiyely compare star formation in the merger to that in the unperturbed galaxies. 
The merging galaxies show a strong starburst, in sharp contrast to their isolated progenitors. Most 
star clusters form in the tidal features. With a mass range of 5 x 10^-5 x lO^M©, they are identified 
as globular clusters. The merger remnant is an elliptical galaxy. Clusters with different mass or age 
haye different radial distributions in the galaxy. Our results show that the high specific frequency 
and bimodal distribution of metallicity obseryed in elliptical galaxies are natural products of gas-rich 
mergers, supporting a merger origin for the ellipticals and their globular cluster systems. 
Subject headings: galaxies: interactions — galaxies: star clusters — galaxies: starbursts — galaxies: 
eyolution — stars: formation 



1. INTRODUCTION 

Interactions and mergers driye galactic eyolution and 
cause starbursts. In the past decade, the Hubble Space 
Telescope (HST) has unyeiled eyidence for the forma- 
tion of massiye star clusters or young gl obular clusters in 
merg ing galaxies, such as NGC 7252 l |Whitm ore et al.l 
1993), the Ante nnae Galaxy ( Whitmore & Schweizcr 



1991. NGC 3597 llHoltzman et al.iil9 96'). and NGC 3921 
ijSchweizer et al.ll996|) . as reyiewed by Whitmore ( 2001). 

There haye been seyeral theoretical models for 
the formation of globular clusters (GCs), includ- 
ing primary models where c lu sters form before 
galaxies (e.g iPeebles fc Dickel Il968j) . secondary 

models wh ere clu sters form with gala xies (e.g. 

iFall fc Reei Il985t iHarris fc Pudrit^ [l^, tertiary 
models whe re clusters form afte r galaxies, such as 
in r nergers (lAshman fc Zep 3 I1992D . and unified mod- 
els l|Elrnegre eiifcE freinovni Models haye been 



reyi ewed by 



lAshman fc Zerl Ip^ . 
and iCarney fc HarrisI l|2001() . The HST obseryations 
haye prompted particular interest in the merger sce- 
nario. There haye been many numerical sim u lations 
of gal axy mergers iToomre fc Toomrd Il972t iWhitf 



iki fc Shapirol Il982t iBarnes ^ Hernouist 
iisiJ'1992': 'Mihos fc Hernquist 1994', 



1^22 jF^^S^ki 

1992^ 'Hernauisi^ '1992'; "Mihos fc Hernquist 
Dubinski et al. 1999: Soringel fc White 1999; Soringel 
200nt iNaa b fc Burkcrt 2001; Barnes 2002), but only 
a few focus on GC formation in mergers. These 
include semi-analytical models ijBeasley et al.l l2002|) . 



combined N-body and smoothed particle hydrodynamics 
(SPH) simulations (,Bekki fc Chiba 2002; Bekk i et all 
l2002t). and adaptiye-gri d cosmological simulations 
ijKraytsoy fc GnedinI 120031) . In these models indiyidual 
GCs are not identified and directly followed. In addition, 
the star formation history of the merger is not compared 
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to that of the isolated progenitors. 

We here model GC formation in both isolated disk 
galaxies and their mergers using N-body/SPH simula- 
tions. Absorbing sink particles are used to directly rep- 
resent indiyidual massiye young star clusters. In § 2 we 
briefly describe our computational method. We discuss 
the burst of GC formation during the merging eyent 
in § 3, and compare with the eyolution of the isolated 
galaxy. We then focus on the cluster mass and age dis- 
tribution at the end of the merger in § 4, and summarize 
in § 5. 

2. NUMERICAL METHOD 

We use the publi cly ayailable SPH code GADGET 
ijSnringel et al.l200l[) , and implement absorbi ng sink par- 
ticles that do not interact hydrodynamically ((Bate et "all 
Il995|) to directly measure the mass of grayitationally 
collapsing gas. To represent GCs, we allow sinks 
to form when the local density exceeds 1000 cm~'^ 
ijBromm fc Clar ke 2002). Our galaxy model consists of 
a dark matter halo, and a disk of star s and isotherma l 
gas. We follow the analytical work b y iMo et alJ lll998l^ 
and the numerical implementation by lSprineel fc Whitd 



((1990) and Soringel (2000). Our simulations meet 
three numerical crit eria, the Jeans resolution crite- 
rion l)Bate fc Burker3fl997|) . the grayit y-hydro balance 
criterion for grayitational softening l|Bate fc BurkertI 

f9^, and the equipartition criterion for particle masses 
tcinmctz & White 1997). We adopt a halo concentra- 
tion parameter C = 5, a spin paramete r A = 0.05, an d 
Hubble constant Hq = 70 km s~^ Mpc'^ ()Springe]|2000|) . 

The galaxy model studied here initially has rotational 
yelocity V200 = 100 km s~^ at the yirial radius where 
the oyerdensity is 200, and a yirial mass of M200 — 3.3 x 
lO^^M©. Total mass in the disk and in the disk gas as 
, , frsjqtio" of the yirial mass are fd = 0.05 and fg = 0.01, 

rklessenOaip^de i « ■ ii i i- r i i '^i i 

rcGpcctiyaly. An isothermal equation 01 state with sound 
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speed Cs = 6 km s^^ is used. The total particle number 
for each single galaxy is iVtot = 1 x 10^, with A^g = 5x10^, 
A'^d = 2 X 10^, and Ni, — i x 10^ for disk gas and stars, 
and halo dark matter, respectively. The corresponding 
gravitational softening lengths are eg — 0.01 kpc, = 
0.1 kpc, and eh — 0.4 kpc. The spatial resolution of 
the gas is thus 10 pc, and the mass per gas particle is 
6.6 X 10"^ M0. This is two orders of magnitude better 
resolution than the models of Bekki et al. (2002). 

We perform two simulations. One follows the evolution 
of the isolated galaxy described above, and the other is 
an equal-mass, head-on merger of two such galaxies. The 
two galaxies are on a parabolic orbit and are separated 
initially by a distance of 300 kpc. In the merger run, 
A'^tot = 2 X 10^, with A'^g = 1 X 10^ gas particles and 
corresponding other particle numbers. The simulations 
are followed for up to 5 Gyr. 

3. STARBURSTS IN THE MERGING EVENT 
3.1. Where Do Globular Clusters Form? 

Our single galaxies begin marginally stable to gravi- 
tational instability (Rafikov 200l|), with initial Toomre 
instability parameter for gas Q — lA. Stars form slowly 
but steadily, mostly in the center. In the merger event, 
the galaxies first collide at t ~ 1.4 Gyr, then separate to 

200 kpc at t ^ 2 Gyr, fall back, and finally merge at 
t ^ 4.5 Gyr. Vigorous starbursts occur in each of the 
two close encounters. 
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Fig. 1. — GC's formed in the merger. The blue- white image 
is the gas surface density at t = 1.53 Gyr, with values given by 
the color bar. The colored dots represent (yellow) older GCs with 
lifetimes Tgc > 10 Myr, (green) young GCs with Tgc < 10 Myr and 
Mgc < 10'' Mq, (red) young GCs with Mgc > 1O''M0, and (black) 
the most massive young GC. 



Most clusters form during the first close encounter at 
t =1.4-1.7 Gyr. To estimate their stellar masses, we 
make the assumption that individual sink particles rep- 
resent regions of dense mole cular gas. Obs ervations by 
iRownd fc YotI^ l|T99ai and iWong &: Blit3 (2002.) sug- 
gest that the local star formation efficiency in molecular 



clouds remains roughly constant. We assume that 35% 
of the mass of each sink particle forms stars, and will 
show in future work that a range of 20-50% is consistent 
with the observed Schmidt law. Using this assumption, 
the cluster mass range is 5 x lO^M© to 5 x 1O^M0. We 
therefore identify these clusters as newly formed GCs. 

Figure ^ shows that during the first encounter GCs 
form in the extended tidal features where the gas 
reaches high density, in agr e ement with obser vations 
(jWhitmore fc Schweize J IT995t iZhang et alJl2001() . The 
derived n iass range agrees we ll with spectroscopic esti- 
mates hv lMenge! et"an l|200l of some young clusters in 
the Antennae Galaxy. 

3.2. Comparison Between Merger and Isolated Calaxy 

Figure|2compares star formation in the isolated galaxy 
with the galaxy merger. In the isolated galaxy, stars 
start to form at t ~ 0.3 Gyr when gravitational insta- 
bility increases the gas density. The cluster number A'gc 
and mass Mgc slowly increase as the galaxy maintains 
marginal stability. During the merging event, the gas 
density is locally increased by shocks and tidal forces, 
leading to starburst behavior. The cluster formation rate 
dMgc/dt increases dramatically, particularly in the first 
encounter when more gas is available. At the end of the 
simulation, A"gc = 73 for the isolated galaxy and 417 for 
the merger. 

The specific frequency ()Harris fc van den BerdflQSll) 
of GCs in a galaxy is Sn = NtlO°-^^^-+'^^\ where A^t is 
the total number of GCs, and the V-band absolute mag- 
nitude of the host galaxy. My = Mq — 2.5\og{Ly/ Lq). 
We use the method of Bekki et al. ( 2002) to estimate 
S'n in our models by assuming M^/Ly 3, as in 
the Milky Way. In order to accoun t for th e clus- 
ter destruction mechanism s des cribed bvlSoitzerl l[l9 871: 
IGnedin fc Ostrikeil llT99l . and 'Fall fc Zhand ll2001f) . we 
assume a 10% survival rate of our clusters. Our final de- 
rived value of S'n depends linearly on this. Each single 
galaxy has a stellar mass of ~ 1.65 x 10^ Mq, so we 
can derive A4v ^ —17.1 mag, so S'n ~ 1-0. Similarly, 
the elliptical merger remnant has Af, ^ 3.3 x IO^Mq, 
which gives Mv ~ —17.9 mag, and Sn ^ 3.0. These 
values are consistent with observations. For example, 
the Milky Way-like spir al M31 has S'n ~ 1.2 ± 0.2 
(jBarmbv fc Huchral 120011), while typi cal ellipticals have 
Sm ~ 3.5 ljHarrislll99H 'McLaughlin' 199^. These es- 
timates do depend on the st ellar mass-to- light ratio 
of th e galaxies (e.g. iLarsen et al.l i2002t |Rhode & Zegfl 
12004!) and survival rate of the clusters ( Soitzer ^983 
Gnedi n fc Ostrikeilll997t IFall fc Z hang 2001; Wh itmOT^ 
i2004^) . Nevertheless, they reflect the difference between 
isolated disk galaxies and mergers, and suggest that 
mergers can explain the considerably higher specific fre- 
quency of GCs observed in elliptical galaxies. 

4. MERGER REMNANT 

Star formation during the merger leads to rapid gas 
depletion. By t = 5 Gyr the two merging galaxies have 
transformed into an elliptical galaxy, whose stellar distri- 
bution is shown in Figure |21 Three distinct populations 
of GCs form at different epochs, those formed during the 
two close encounters are almost coeval, as they form over 
a very short period during the encounters. 
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Fig. 2. — Number A^gc (top) and fractional mass A'/gc / -A/gas 
{bottom) of GCs over time compared between two merging galaxies 
{solid lines) and the same two galaxies isolated {dashed lines). Note 
the sudden jumps in the merging galaxies, representing starbursts. 



FigureOlshows the GC system is centered on the steUar 
spheroid. Massive GCs are concentrated towards the cen- 
ter of the eUiptical galaxy that remains after the merger, 
while less massive ones are found at greater radii as well. 
Old GCs formed prior to the first encounter and the 
youngest GCs formed in the second encounter are con- 
centrated towards the center of the galaxy, while GCs 
formed in the first encounter are also found at large radii. 
The oldest GCs in our simulation are also the most mas- 
sive, and quickly sink towards the center by dynamical 
friction. The youngest ones occur in the galactic center 
because remaining gas assembles there during the final 
merger. During the first encounter, on the other hand, 
GCs form in the tidal tails with greater angular momen- 




FlG. 3. — Log of stellar surface density in the merger remnant, 
with values shown in color bar, compared to distribution of GCs 
{dots) in mass {top panel) and age {bottom panel). Mass bins shown 
are Afgc < lO^M© {blue), 10^ < Mgc < {cyan), lO'^-^ < 

Mgc < lO'^M© {green), Mgc > lO'^M© {red), and the most massive 
GC, Mgc = 7.8 X IO^Mq {black). Age bins shown are GCs formed 
at t < 1.4 Gyr during initial approach before the merging event 
begins {blue), GCs formed at 1.4 Gyr < t < 1.9 Gyr during the 
first encounter {green), and GCs formed at 4.0 Gyr < 4 < 4.6 Gyr 
during the second encounter {red). 



tum, and so have a broad radial distribution. 

If we assume a metallicity-ag e relation such that 
older GCs have lower metallicities l)Rich et al.ll2nniD . we 
clearly see a bimodal metallicity distribution. Young, 
high-metallicity GCs are centrally concentrated, while 
older, lower-mctallicity GCs ex tend to larger radii, 
in ag reement with observations l)Diorgovski fc MevlanI 
Il994|) . The isolated galaxy forms clusters slowly and 
steadily, so age, as well as metallicity, follow a smooth 
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distribution. There is no bimodality. Our r esults thus 
support the arguments of lAshman fc ZepJ l)1992H and 
iKundu fc Whitmord l)200H) that mergers are required to 
produce a bimodal distribution. 

5. SUMMARY 

We present high-resolution simulations of star cluster 
formation in both a single galaxy and a major merger, 
using a three-dimensional SPH code that includes ab- 
sorbing sink particles to represent massive star clusters. 
This allows direct identification of individual clusters and 
tracking of their orbital evolution over several gigayears. 
The merging galaxies show bursts of massive star cluster 
formation, in sharp contrast to the steady but slow for- 
mation in an isolated galaxy. Most new clusters form in 
the tidal tails and bridges between the merging galaxies. 
They are identified as progenitors of globular clusters, 
although it should be emphasized that we do not include 
dynamical destruct ion of clusters (e.g. iSpitzeil Il987t 
iCxnedin fc Ostrikcr 1997": "Fall fc Zhang 20011 in the sim- 
ulations. Dynamical destruction will change the number 
and mass of GCs, affecting our derived specific frequency 
(where we accounted for it very crudely) and spatial dis- 
tribution. Inclusion of a destruction model to determine 
the evolution of our cluster particles will be necessary to 
fully address these questions. 

The estimated specific globular cluster frequency S'n 
in the elliptical galaxy resulting from the merger exceeds 
by a factor of 3 that in an isolated galaxy with the same 
parameters as the merging galaxies, most of the enhance- 
ment being in metal-rich GCs. This supports the idea 



that the higher S'n of metal-rich GCs observed in ellipti- 
cals is produced by mergers. However, elliptical galaxies 
also show higher S'n for metal-poor GCs. Early mergers 
could explain t his, though other me chanisms have also 
been proposed (Rhode & Zcpf '20(i^. Clusters formed 
during different phases of the merger have distinct radial 
distributions. The spatial distribution of metal-poor old 
clusters formed before the first encounter as well as of 
metal-rich younger clusters from the final encounter is 
centrally concentrated, while clusters formed in the first 
encounter are more widely dispersed. This suggests the 
observed multimodal metallicity distribution of globular 
clusters in elliptical galaxies is the direct result of merger 
processes. 
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